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We discuss duality in "two-photon" -like processes in the scalar <p% model and also in the process 
7*7 — > 7T7T in QCD. Duality implies the equivalence between two distinct nonperturbative mech- 
anisms. These two mechanisms, one involving a twist-3 Generalized Distribution Amplitude, the 
other employing a leading-twist Transition Distribution Amplitude, are associated with different 
regimes of factorization. In the kinematical region, where the two mechanisms overlap, duality is 
observed for the scalar ip% model, while in the QCD case the appearance of duality turns out to be 

in : 

sensitive to the particular nonperturbative model applied and can, therefore, be used as a tool for 
selecting the most appropriate one. 
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INTRODUCTION 



The only known method today to apply QCD in a rigorous way is based on the factorization of the dynam- 
ics and the isolation of a short-distance part that becomes this way accessible to perturbative techniques of 
quantum field theory (see, [J 0, Q and for a review, for instance, @| and references cited therein). Then, the 
' conventional systematic way of dealing with the long-distance part is to parameterize it in terms of matrix 
(~~. , elements of quark and gluon operators between hadronic states (or the vacuum) . These matrix elements stem 
lO ' from nonperturbative phenomena and have to be either extracted from experiment or be determined on the 
lattice. In many phenomenological applications they are usually modeled in terms of various nonperturbative 
■ methods or models. 

Generically, the application of QCD to hadronic processes involves the consideration of hard parton sub- 
CT\ . processes and (unknown) nonperturbative functions to describe binding effects. Prominent examples are hard 
exclusive hadronic processes which involve hadron distribution amplitudes (DAs), generalized distribution am- 
plitudes (GDAs), and generalized parton distributions (GPDs) 0,0, HQ- Applying such a framework, collisions 
of a real and a highly- virtual photon provide a useful tool for studying a variety of fundamental aspects of QCD. 

Recently, nonperturbative quantities of a new kind were introduced — transition distribution amplitudes 
(TDAs) [3, [HI [H| — which are closely related to the GPDs. In contrast to the GDAs, the TDAs appear in 
the factorization procedure when the Mandelstam variable s is of the same order of magnitude as the large 
photon virtuality Q 2 , while t is rather small. Remarkably, there exists a reaction where both amplitude types, 
GDAs and TDAs, can overlap. This can happen in the fusion of a real and transversely polarized photon with a 
highly-virtual longitudinally polarized photon, giving rise to a final state which comprises a pair of pions. The 
key feature of this reaction is that it can potentially follow either path: proceed via twist-3 GDAs, or go through 
the leading-twist TDAs, as illustrated in Fig. [TJ Such an antagonism of alternative factorization mechanisms in 
this reaction seems extremely interesting both theoretically and phenomenologically and deserves to be studied 
in detail. 

The intimate relation between these two mechanisms in the production of a vector-meson pair was analyzed 
in [l2| and it was found that these mechanisms can be selected by means of the different polarizations of the 
initial-state photon. In contrast, for (pseudo)scalar particles, such as the pions, this effect is absent enabling us 
to access the overlap region of both mechanisms and their duality as opposed to their additivity. 

In this talk, we will report on the possibility for duality between these antagonistic mechanisms of factoriza- 
tion, associated either with GDAs or with TDAs, in the regime where both Mandelstam variables s and t are 
rather small compared to the large photon virtuality Q 2 . 
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Figure 1: Two ways of factorization: via the GDA mechanism and via the TDA mechanism. 



II. REGIMES OF FACTORIZATION WITHIN THE ^-MODEL 

Consider first the factorization of the scalar ip\ model in Euclidean space. To study the four-particle amplitude 
in detail, it is particularly useful to employ the a-representation — see Q ■ Then, the contribution of the leading 
"box" diagram can be written as (while details can be found in [13]) 



A(s, t, to 2 ) = 



16tt 2 
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n dcti 

D 2 



■ exp 
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sa2dA + toi«3 + m 2 D 2 ) 
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where m 2 serves as a infrared (IR) regulator, s > 0, t > are the Mandelstam variables in the Euclidean 

4 

region, and D — ^ o^. Assuming that q 2 = Q 2 is large compared to the mass scale to 2 (which simulates here 

the typical scale of soft interactions), the amplitude (p} can indeed be factorized. As regards the other two 
kinematic variables s and t, one can identify three distinct regimes of factorization: (a) s <C Q 2 while t is of 
order Q 2 ; (b) i < Q 2 while s is of order Q 2 ; (c) s, t <C Q 2 . 

Regime (a): The process is going through the s-channel. In this regime, the main contribution in the integral 
in Eq. (|TJ) arises from the integration over u\ when ot\ ~ 0: 



A^ DA {s,t,m 2 ) 



16tt 2 
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(2) 



Schematically this means that the propagator, parameterized by a%, can be associated with the partonic (hard) 
subprocesses, while the remaining propagator constitutes the soft part of the considered amplitude, i.e., the 
scalar version of the GDA. 

Regime (b): Here we have to eliminate from the exponential in Eq. ([I]) the variables Q 2 and s, which are 
large. This can be achieved by integrating over the region a-i ~ 0. Performing similar manipulations as in 
regime (a), we find that the scalar TDA amplitude can be related to the scalar GDA via _4^ s DA (s, t, m 2 ) = 
A^ BA (t,s,m 2 ). 

Regime (c): The relevant regime to investigate duality is when it happens that both variables s and t are 
simultaneously small compared to Q 2 , i.e., when s, t <C Q 2 ■ In this case, there are two possibilities to extract 
the leading Q 2 -asymptotics, notably, we can either integrate over the region a\ ~ 0, or integrate instead over 
the region ~ 0. Clearly, these two options can be associated with (i) the GDA mechanism of factorization 
with the meson pair scattered at a small angle in its center-of-mass system or, alternatively, (ii) with the 
TDA mechanism of factorization. We stress that we may face double counting when naively adding these 
two contributions. We interpret such a behavior as a signal of an ingrained tendency for duality between the 
GDA(s-channel) and the TDA (t-channel) factorization mechanisms. 

In order to verify the appearance of duality we carry out a numerical investigation of the exact and the 
asymptotic amplitudes. In doing so, we introduce the following ratios R\ — _4^ S DA /.A and i?2 = Aq^/A. 
Appealing to the symmetry of these ratios under the exchange of the variables s <-> t, we take t/Q 2 to be 0.01 
and look for the variation of the ratios with s/Q 2 . This variation is illustrated in Fig.^from which one sees that 
in the region where s/Q 2 is rather small, i.e., in the range (0.01, 0.05), both asymptotic formulae are describing 
the exact amplitude with an accuracy of more than 90%. This behavior supports the conclusion that, when 
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Figure 2: The ratios Ri and R2 as functions of s/Q 2 



both Mandelstam variables s/Q 2 and t/Q 2 assume values in the wide interval (0.001, 0.7), duality between the 
TDA and the GDA factorization mechanisms emerges. 



III. TDA- AND GDA- FACTORIZATIONS FOR 77* tvtv 

Having discussed the appearance of duality between the GDA and the TDA factorization schemes within a toy 
model, we now turn attention to real QCD. To analyze duality, we consider the exclusive Tr + n~ production in a 
7t7l collision, where the virtual photon with a large virtuality Q 2 is longitudinally polarized, whereas the other 
one is quasi real and transversely polarized. Notice that the GDA and the TDA regimes correspond to the same 
helicity amplitudes. Given that the considered process involves a longitudinally and a transversally polarized 
photon, we are actually dealing with twist-3 GDAs [l4j. On the other hand, for the twist-2 contribution, related 
to the meson DA, we use the standard parametrization of the 7r + -to-vacuum matrix element which involves a 
bilocal axial- vector quark operator Finally, the 7 — * 7r~ axial- vector matrix elements can be parameterized 
in the form, cf. fl"ot | . 

(ir-(p 2 )\TP(-z/2) lal5 [-z/2;z/2}Tp(z/2)\j(q',e')) £ -^-e' T ■ AtP u Ai(x, £, t) , (3) 

Jit 

where P = (p2 + 1 )/2, and A = P2 — q' , and noticing that the symbol = means Fourier transformation and that 
the vector matrix element does not contribute here. To normalize the axial- vector TDA, A%, we express it in 
terms of the axial- vector form factor measured in the weak decay ir — > lv\^j (l3l. [Tol HH]. The helicity amplitude 
associated with the TDA mechanism reads 

c-'(j') AT 

,TDA _ ^-TDA £ (A) 



with 



where 
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(5) 
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TDA^j = J dxAi(x,£, 



"l^-tT=l- 



employing the 1-loop a s (Q 2 ) in the MS-scheme with A QC d = 0.312 GeV for Nf = 3 [Note that there is 

only a mild dependence on Aqcd-] 

Turning now to the helicity amplitude, which includes the twist-3 GDA, we anticipate that it can be written 
as (see, for example, [1J|) 

aGDK _ -T-GDA £ ' j± f 7 \ 



with 

-GDA_ ^1+Q 2 



jrGDA _ 2 + eg) (tw- 3 GDA WW), (8) 



where 



tw 



l 

-3GDAWwj= J dyd^iyX.W 2 )^ " ^) . (9) 
o 



with the partial derivative being defined by 9^ = d/d(2( — 1). In deriving ((8|), we have used for the twist-3 
contribution the Wandzura-Wilczek approximation. Duality between expressions ([5]) and ([5]) may occur in 
that regime, where both variables s and t are simultaneously much smaller in comparison to the large photon 
virtuality Q 2 . More insight into the relative weight of the amplitudes with TDA or GDA contributions can be 
gained once we have modeled these non-perturbative quantities. We commence our analysis with the TDAs 
and, assuming a factorizing ansatz for the i-dependence of the TDAs, we write A\{x, £, t) — 2-^- FA{t)Ai(x, £), 
where the ^-independent function Ai(x,^) is normalized to unity. To satisfy the unity-normalization condition, 
we introduce a TDA defined by 

Al{Xjl) = _^i £i!L_ (io) 

/ dxA™*-™™^, 1) 
-l 

and continue with the discussion of the i-independent TDAs. Recalling that we are mainly interested in TDAs 
in the region £ = 1 [l], [2j , it is useful to adopt the following parametrization 



A™ n -™ Im (x, 1) = (1 - x z ) \\ + ai C[ J, ^{x) + a 2 C y 2 °<'>{x) + a 4 Q 0/ >)J , (11) 

where a\, a,2, 04 are free adjustable parameters, encoding nonperturbative input, and the standard notations for 
Gegenbauer polynomials are used. It is not difficult to show that the TDA expressed by Eq. (jlip results from 
summing a D-term, i.e., the term with the coefficient a\, and meson-DA-like contributions. For our analysis, 
we suppose that a\ = do 0j which is equal to —0.5 in lattice simulations. With respect to the parameters 0,2 
and (Z4, we allow them to vary in quite broad intervals, notably, 0,2 G [0.3, 0.6] and 0,4 £ [0.4, 0.8], that would 
cover vector-meson DAs with very different profiles at a normalization scaleu 2 ~ 1 GeV 2 (see, for example, 
[l8|). The function $i(z,£) is rather standard and well-known (details in I5l [l3l '). 

We close this section by summarizing our numerical analysis presented in [131 ] . We calculated both functions 
jrTDA anc j ^rGDA^ anc j s ] 10w results in Fig. [3l The dashed line corresponds to the function JF TDA , where we 
have adjusted the free parameters to 0,2 — 0.6, — 0.8. The results, obtained for rather small values of these 
parameters, are displayed by the broken lines in the same figure. The dotted line denotes the function 
with 02 = 0.5 and 04 = 0.6, whereas the dashed-dotted line employs a 2 — 0.3 and 0,4 = 0.4. For comparison, we 
also include the results for ^ rGDA . In that latter case, the dense-dotted line corresponds to the GDA amplitude, 
where the expression for B\2 bas been estimated via Eq. (20) of [l3j], while the solid line represents the simplest 
ansatz for B\ 2 with = 0.5. From this figure one may infer that when the parameter -B12, which parameterizes 
the GDA contribution, is estimated with the aid of the Breit-Wigner formula (provided s, t <§; Q 2 ), there is 
duality between the GDA and the TDA factorization mechanisms. Hence, the model for $i(z,£), which takes 
into account the corresponding resonances, can be selected by duality. 
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Figure 3: Helicity amplitudes J rTUA and J^ GUA as functions of Q 2 
value of s/Q 2 varies in the interval [0.06, 0.3]. 



using a\ = —0.5 found in lattice simulations. The 



IV. CONCLUSIONS 

We have provided evidence that when both Mandelstam variables s and t turn out to be much less than 
the large momentum scale Q 2 , with the variables s/Q 2 and t/Q 2 varying in the interval [0.001, 0.7], the TDA 
and the GDA factorization mechanisms are equivalent to each other and operate in parallel. We have also 
demonstrated that duality may serve as a tool for selecting suitable models for the nonperturbative ingredients 
of various exclusive amplitudes entering QCD factorization. In this context, we observed that twist-3 GDAs 
appear to be dual to the convolutions of leading- twist TDAs and DAs, multiplied by a QCD effective coupling. 
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